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Abstract 


Dwarf irregular galaxies (dIrrs), as rotationally supported systems, have more reliable J-factor measurements than 
dwarf spheroidal galaxies and have received attention as targets for dark matter detection in recent years. In this 
paper, we use 10 yr of IceCube muon-track data and an unbinned maximum-likelihood-ratio method to search for 
neutrino signals beyond the background from the directions of seven dlIrrs, aiming to detect neutrinos produced by 
heavy annihilation dark matter. We do not detect any significant signal. Based on such null results, we calculate the 
upper limits on the velocity-averaged annihilation cross section for 1 TeV-10 PeV dark matter. Our limits, 
although weaker than the strictest constraints in the literature in this mass range, are also a good complement to the 
existing results considering the more reliable J-factor measurements of dlrrs. 
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1. Introduction 


Astrophysical and cosmological evidence has suggested the 
existence of particles beyond the Standard Model, dark matter 
(DM), which account for about 84% of the matter component 
of the universe (Planck Collaboration et al. 2016, 2020; 
Bertone & Hooper 2018; Arbey & Mahmoudi 2021). Although 
its nature is still unknown, many models, such as primordial 
black holes (PBH) (Carr & Kuhnel 2021), weakly interacting 
massive particles (WIMPs) (Bertone et al. 2005; Roszkowski 
et al. 2018), axion/axion-like particles (ALPs) (Saikawa 2017) 
and sterile neutrinos (Kopp 2021) are proposed as dark matter 
candidates. In recent years, most dark matter searches have 
focused on the mass range of the sub-eV (e.g., ALPs) or the 
GeV -TeV scale (e.g., WIMPs) and have not found a significant 
signal (Acharyya et al. 2023). The situation prompts us to 
consider the models in a wider mass range, one example is the 
ultra-heavy dark matter (UHDM). 

The UHDM has a mass of above 10 TeV, which means it is 
too massive to be produced in current colliders. There are 
already many viable candidates for UHDM (Chung et al. 1998; 
Berlin et al. 2016a, 2016b; Kolb & Long 2017; Carney et al. 
2022), and in this paper we focus on the annihilation UHDM 
which has received less attention before. In the scenario of 
thermal production of dark matter, the cross section of DM 
annihilation to Standard Model particles scales as My ? (Tak 
et al. 2022). As the mass increases, the cross section generally 
decreases. If it becomes too small, then the DM will not be 
sufficiently depleted before freeze-out, leading to an incon- 
sistency with the observed cosmological DM density. This 
unitarity consideration limits the DM mass to be below 


7100 TeV (Griest & Kamionkowski 1990; Bramante & 
Unwin 2017). Although this unitarity boundary has long 
hindered people's interest and enthusiasm in searching for the 
annihilation UHDM, several mechanisms have been proposed 
that can violate this boundary (Kim & Kuflik 2019; Carney 
et al. 2022). For instance, if the DM exists in a composite state 
rather than a point-like particle, thermal DM with PeV-scale 
masses could be possible (Harigaya et al. 2016; Tak et al. 
2022). In such a scenario, the unitarity bound can be violated 
because the DM annihilation cross section no longer scales as 
My 2 Instead, it is now related to the geometric size of the 
composite DM. 

One way to detect dark matter particles is to observe the 
gamma-rays, neutrinos from the DM annihilation or decay in 
the astrophysical environments. Promising sites for the 
observation of neutrinos from dark matter annihilation include 
the galaxy clusters, Galactic halo and center, as well as the 
cores of the Sun and Earth (Zornoza 2016). Among other 
targets, dwarf spheroidal galaxies (dSphs), one type of the 
Milky Way satellite galaxies, are also considered as promising 
objects because of their relatively close proximity and clean 
astrophysical background (Wood et al. 2015). However, the 
stellar kinematics of these pressure-supported galaxies are 
highly uncertain, which makes it difficult to determine their 
DM profiles very well (Battaglia et al. 2013). The accurate 
determination of the DM profile is very important for 
estimating astrophysical factor (J-factor) of DM indirect 
detection (Hayashi et al. 2016). 

Dwarf irregular galaxies (dIrrs) are targets that we have paid less 
attention to Hernández Cadena et al. (2019), Abdallah et al. (2021), 
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Gammaldi et al. (2021), Alfaro et al. (2023), and they have some 
properties similar to dSphs, such as high DM abundance and low 
astrophysical background. Unlike dSphs, dIrrs are located further 
away and they are rotationally supported star-forming dwarf 
galaxies (Oh et al. 2015; Gammaldi et al. 2021) (in contrast, dSphs 
are pressure-supported systems). Their DM profiles can be 
obtained from their rotation curves. Meanwhile, dlrs are in 
general larger galaxies compared to dSphs and have larger dark 
matter halos (their halo masses are M200 = 10’-10'°M.), and may 
host more DM subhalos (Karukes & Salucci 2017). The standard 
ACDM cosmological model predicts that dark matter halos are 
filled with many self-bound substructures, which has been 
confirmed by N-body simulations (Frenk & White 2012; Planck 
Collaboration et al. 2016). The presence of subhalos is expected to 
significantly enhance the DM signal (Lacroix et al. 2022; Runburg 
et al. 2022). Therefore, the enhancement factor due to subhalos 
increases the detectability of the dIrrs. These relevant properties 
make dlrrs also interesting targets for DM searches. 

Previous works by IceCube collaboration, which have searched 
for dark matter signals from the Galactic center (Albert et al. 
2020b; Iovine & IceCube Collaboration 2021), Galactic halo 
(Abbasi et al. 2023; Desai et al. 2023), Sun center (Argüelles et al. 
2019; Abbasi et al. 2022), Earth center (Renzi 2019; Saveliev & 
Hyde 2022), and nearby galaxies and galaxy clusters (Aartsen 
et al. 2013a; Jeong et al. 2022), found no significant excess 
beyond background expectations. In our previous work, we have 
used IceCube observations on dwarf spheroidal galaxies to search 
for dark matter (Guo et al. 2023). In this paper, we further search 
for DM in dwarf irregular galaxies on the basis of the previous 
work. Different from the previous analysis, we analyze dwarf 
irregular galaxies as extended sources and consider the enhance- 
ment of the annihilation signal due to their internal subhalos. We 
will search for UHDM induced neutrinos from dIrrs using 10 yr 
(2008 April to 2018 July) of Icecube muon track data (IceCube 
Collaboration et al. 2021) and an unbinned maximum-likelihood- 
ratio method (Braun et al. 2008). We select seven dlrrs, four 
sources (WLM, NGC 6822, IC10, IC1613) are the most 
promising objects for DM detection as proposed in previous 
works (Gammaldi et al. 2018, 2021), and the remaining three 
sources are selected from (Alfaro et al. 2023). Note that the 
IceCube 10 yr muon track data have a much higher sensitivity to 
the northern hemisphere sky. Nevertheless, we have also included 
two sources in the southern sky (WLM and NGC 6822) in the 
sample to demonstrate this sensitivity difference. The sample 
information is shown in Table 1. 


2. Signal Expectations 
The expected differential neutrino flux for the annihilation of 
dark matter particles is calculated as follows: 
dọ, ov) dN, 
(E) = ) 
dE, 8amy dE, 


(E Jao x B (1) 
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where (ov) is the velocity-averaged annihilation cross section, 
my is the mass of the dark matter particle, dN,/dE, is the 
annihilation spectrum of DM particles, and B is the boost 
factor. 

We assume dark matter annihilates into the standard model 
particles with a 10096 branch ratio and calculate the energy 
spectrum of neutrinos dN, /dE,. For this calculation, we use the 
HDMSpectra (Bauer et al. 2021), a Python package providing 
spectra of various annihilation and decay channels of heavy 
dark matter. In this work, we consider the DM masses ranging 
from 1 TeV to 10 PeV and six annihilation channels: bb, if, 
T'r,wU,,u' p and WW. Figure 1 shows the neutrino 
spectra for the annihilation of a 100 TeV DM through different 
channels. 

The flux is proportional to the J-factor (the integral over the 
square of the dark matter density), 


nme J ; J __ dl ping") (2) 


where / is the distance along the line of sight (lo.s.) and 
r(I) = 4D? — 21Dcos0 + I? with D the distance to the 
source and 0 the angular separation between the directions of 
the source center and l.o.s. In our work, we adopt the J-factors 
reported in (Alfaro et al. 2023). 

The universal rotation curve model (Persic et al. 1996) can fit 
the dlrr rotation curve well (Gammaldi et al. 2018), and the 
DM distribution in these galaxies is well described by a Burkert 
profile (Burkert 1995): 


3 
a (3) 


Pom) = POE EET 


where po is the normalization density, and rọ is the core radius. 

The substructures within dIrrs can increase the expected DM 
flux. We parameterize this effect with the boost factor B in 
Equation (1). For the calculation of the boost factors we mainly 
refer to the formulas and results in (Moliné et al. 2017). They 
considered the VL-II (Diemand et al. 2008) and ELVIS 
(Garrison-Kimmel et al. 2014, 2017) N-body cosmological 
simulations to study the properties of subhalos, obtaining the 
concentrations of the subhalos as a function of the distance Rub 
to the main halo center. Together with the concentration-mass 
relation derived in (Pieri et al. 2011) for the field halo, we can 
calculate the luminosity Lsmooh(M) of the smooth halo for a 
given halo mass M and the luminosity Lsmooth(™M, Xsub) of the 
subhalo for a given halo mass m at the distance RA, 
(Xsub = Reub/Rvir with Ry, the virial radius). Then the boost 
factor can be calculated with: 


B(M) = 


3 f dN (m) 
L smooth(M) Min dm 


1 
x dm f dxsiol1 + Bom|L(m, xan — 
0 
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Table 1 
Sample of dIrr Galaxies 
Name R.A. (deg) Decl. (deg) D (Mpc) Ovir (deg) os (deg) logy) J (GeV? cm?) logio Moo (Ms) Boost Factor TS 
IC10 5.10 59.29 0.79 3.857 0.078 17.857 10.6 5.53 0.04 
IC1613 16.19 2.13 0.76 2.361 0.120 17.632 10.0 5.01 0.03 
NGC 6822 296.23 —14.80 0.52 5.325 0.159 18.173 10.5 5.45 0.00 
WLM 0.49 —15.46 0.98 2.609 0.082 18.062 9.6 4.72 0.00 
DDO133 188.22 31.54 4.88 0.784 0.033 17.501 10.5 5.45 0.13 
DDO154 193.52 27.15 4.04 0.695 0.030 17.800 10.3 2:25 0.00 
DDO168 198.61 45.91 4.25 1.142 0.032 17.365 10.9 5.86 0.46 


Note. Columns: (1) galaxy name (2) R.A. (3) decl. (4) galaxy distance (5) the angle to virial radii (6) source extension (7) the astrophysical factor for annihilation 


(8) halo mass (9) Boost Factor (10) TS value. 
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Figure 1. Spectra of neutrinos from the annihilation of DM particles with a 
mass of 100 TeV for six channels, assuming that the branching ratio for each 
channel is 100%. 


The dN(m)/dm is the subhalo mass function, which reads 
dN (m)/dm = A/M(m/M) *. The normalization factor is 
equal to A — 0.012 for a slope of the subhalo mass function 
of a=2 and to A=0.03 for a= 1.9 (Sánchez-Conde & 
Prada 2014). It is important to note that there are large 
uncertainties in the calculation of the boost factor B(M), which 
depends on the index a of the mass function, the assumed 
minimum subhalo mass Min, the concentration-mass relation 
adopted for subhalos and whether or not the tidal stripping 
effect being taken into account. In this work, we adopt the 
median value of the boost factors reported in (Moliné et al. 
2017) the mass function index is assumed to be 1.9, the 
minimum halo mass is 10 °M., and the tidal stripping effect is 
not considered. The boost factors are shown in Table 1. 


3. Icecube Data Analysis 


IceCube is a cubic-kilometer neutrino detector located at the 
South Pole. It consists of 5160 optical photosenors, which are 


connected to the data acquisition system via cables ("strings"). 
Neutrinos interact with nuclei and electrons in the ice 
producing charged particles and emitting Cherenkov light, 
which is detected by the photosensors (Aartsen et al. 2017b). 
The geometry and sensitivity of the photosensors lead to an 
effective low-energy threshold of about 100 GeV for neutrinos 
(Yoshida & IceCube Collaboration 2003). 

In this work, we utilize the IceCube updated public data of 
muon tracks. This data set consists of muon tracks observed by 
IceCube from April 2008 to July 2018, including a total of 
1134450 muon-track events. We use all IC40 to IC86-VII data, 
with the number in the name corresponding to the number of 
installed detector strings. The data include direction and energy 
information of neutrino events, as well as instrument response 
functions (IRFs, including effective area, smearing matrix, etc) 
for different periods over the ten years (IceCube Collaboration 
et al. 2021). We process the 19 x 2 double-counted tracks in 
the data set found in (Zhou & Beacom 2022), following the 
procedure in (Chang et al. 2022).! 

The analysis is performed by using an unbinned maximum 
likelihood ratio method (Braun et al. 2008; Abbasi et al. 2011; 
Aartsen et al. 2013b, 2017a; Li et al. 2022; Pasumarti & 
Desai 2022). The likelihood function is defined as follows: 


Lin) i ns gk ( x] k (5) 
Ns) = s i I 5 B; > 
k i=1 N* N* 


where $; and B; are the signal and background probability 
density functions (PDFs), respectively, N * is the total number 
of events in the sample k and nč is the number of signal events 
which is dependent on the model parameters. The superscript k 
refers to the kth IceCube data sample (IC40-IC86VII). 

The signal and background PDFs for the ith event are 
functions of the reconstructed direction x; and the reconstructed 
muon energy £;. For time-integrated searches the signal PDF S; 


! See https:/ /github.com/beizhouphys/IceCube data 2008-2018 double 


counting corrected for the updated data set. 
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is given by: 
Si = se (xi, ajx;) x SUEY): (6) 


where S;?*" is the spatial contribution which depends on the 
angular uncertainty of the event c; and the angular separation 
between the directions of the event x; and the source x,. We 
model the spatial term of the signal PDF as a two-dimensional 


Gaussian, 
l |x; am xj 
exp| - ———— (7) 
no? i | 202 


t 


s" (xi, aixs) = 


where x; is the direction of the ith event, and x, is the direction 
to the source. Given the larger size of the dIrrs (compared to 
dSphs), it is better to analyze them as extended sources to get 
more realistic limits. In Table 1, we present the extension 
corresponding to the virial radius and the 68% containment 
angle of these dIrrs. When searching for an extended source, 
the value of c; is replaced with o = 4/0? + c? (Aartsen et al. 
2014; Albert et al. 2020a), where o, quantifies the extension of 
the source. In this work, we adopt the angle containing 6896 of 
the total emission (Ko = 0.68 x Jot) as c, (Geringer-Sameth 
et al. 2015; Di Mauro et al. 2022). The signal PDF becomes, 


1 lx; — eA 
ex i (8) 
2n(o; + 0j) ( 2(oi + 5) 


spat __ 
5 i 


The energy term of the signal PDF, 57", is a function of the 
reconstructed muon energy proxy £;. Assuming a power law 
spectrum of spectral index y and for a given decl. 6, of the 
source,” the PDF is (Huang & Li 2022), 


SE, Al CE, &) MEE, &)dE, 
[E^ As (EL, ô )dE, 


ener __ 
5; 


(9) 


For the search for DM signals, the E; ? should be replaced by a 
DM spectrum (i.e., Equation (1)). 

The background PDF B; is obtained directly from the 
experimental data and is given by 


BÉ = BEP] (6) x Be" (Eil) (10) 


spat 


The spatial term, B;", is the normalized event number per unit 


solid angle and is varied with the decl. ó (Zhou et al. 2021), 


NÈ; 
N x Asa 


B" (5) = (11) 


where N5,+3 is the number of events in the decl. band of 6; + 3 
and A(J5,3 is its corresponding solid angle, N, is the total 
number of events in the sample k. The background PDF is R.A. 


: Considering the small extension of the dlrrs we are studying (see c; in 
Table 1), we disregard the variations of the IRFs within the region covered by 
the source, and adopt the instrument response toward the center of the source 
(as, ôs). 
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independent because the IceCube observatory is located at the 
South pole. 
The B” represents the energy term of the background PDF, 
k 


NE 
BS" (Eê) = ——— (12) 
N5,43 AE 


where Nj is the number of events within the decl. range of 
6 € [6; — 3°, 6; + 3°) and the muon energy range of 
logo Erec € [log Ej, log E; + 0.1) for the sample k. 

We maximize the likelihood in Equation (5) to derive the 
best-fit model parameters. A likelihood ratio test statistic of the 
signal is defined as, 


TS — 21oe| Lf) | (13) 


L(n, = 0) 


where n, is the number of signal events, and 7, is the best-fit 
value of the parameter. 

For searching for DM signals, the number of signal neutrinos 
from DM annihilation in Equation (5) is related to the DM 
parameters through, 


Emax 
m= Tive f Aa DEE a 
E, 


min dE, 
where Tiive is the detector uptime, A, is the effective area of 
the detector, d, /dE,(E,) is the energy spectrum of dark 
matter annihilating into neutrinos given by Equation (1), and 
the Emin and Emax are set to the minimum and maximum 
energies allowed by the effective area provided in the 10 yr 
IceCube data. 


4. Result 


Assuming a power-law spectrum, the TS values of the 
putative neutrino signals from these dlrrs are listed in Table 1. 
We do not find any significant excess beyond the background 
expectation in the directions of the seven sources within our 
sample. So we derive the flux upper limits at the 95% 
confidence level, assuming a source energy spectrum of E ?, 
E ?? and E 5, respectively. The results are shown in Figure 2. 

Next, we explore the neutrino signals from the annihilation 
of dark matter. We use Equation (1) to perform searches for the 
DM mass range from 1 TeV to 10 PeV and different 
annihilation channels, and find no annihilation signal. We 
calculate the upper limits at the 95% confidence level in the 
(ov) — m, parameter space for each dIrr. In Figure 3 we show 
the limits for the six annihilation channels (bb , ff , TT, Va Vs 
p's, W'W ). The strongest constraints that we obtain are for 
the xx — wtp” channel. 

In Figure 3, for the two southern sky sources, it can be seen 
that they place overall much weaker limits on the (ov) than the 
dIrrs in the northern sky. However, at large DM masses 
(> 10° GeV), they give stronger constraints. Such results are 
reasonable. In the southern hemisphere, since the event rate of 
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Figure 2. The 9596 C.L. upper limits on the neutrino flux from the seven dlrrs 
at 100 TeV assuming a source spectrum of E". 


atmospheric muons is orders of magnitude higher than that of 
neutrinos (Aartsen et al. 2017b), harsh cuts on the reconstruc- 
tion quality and minimum energy are applied in order to 
suppress this atmospheric muon background, and the cuts 
remove most southern hemisphere events with an estimated 
energy below ~10 TeV (IceCube Collaboration et al. 2021) (so 
for NGC 6822 and WLM we only show limits for DM masses 
of >10 TeV in Figure 3). In contrast, for the northern 
hemisphere, we can filter out the atmospheric muon contam- 
ination by the Earth. For neutrinos at the high-energy end 
(>10° GeV), on the other hand, events from the northern sky 
may interact far from the detector while crossing the Earth, 
resulting in unobservable energy losses and thus difficult to 
reconstruct. This leads to a significantly smaller effective area 
for the northern sky than the southern sky at 10$ GeV 
energies (one can see Figure 2 of IceCube Collaboration et al. 
2021). Since the energy spectra of the atmospheric muon and 
neutrino fluxes are expected to be softer than the astrophysical/ 
DM neutrinos, their contribution to the background becomes 
less significant at these high energies. Therefore, for the large 
mass DM, the analysis of the two dlrrs (NGC 6822, WLM) 
gives more strong constraints on the (cv). 

To improve the sensitivity, we also jointly analyze the data 
from the five northern-sky sources simultaneously. In this case, 
the combined likelihood function is 


É((ov), my) = I] £ j(ns((ov), My, Jj)) (15) 
J 


with £;(n;) the likelihood in Equation (5) for the jth source. 
The details of the method of the combined analysis can be 
found in (Guo et al. 2023; Zhu et al. 2023). The limits from the 
combined analysis are shown as black lines in Figure 3. The 
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results of the combined analysis are about 2.5 times stronger 
than the ones by DDOI54, which gives the strongest 
constraints among the seven dlir galaxies. 

Finally, in Figure 4 we compare the results of this analysis 
with the limits obtained in other works (Aartsen et al. 2013a; 
Ackermann et al. 2015; Abdallah et al. 2021; Gammaldi et al. 
2021; Miener et al. 2022; Acharyya et al. 2023; Alfaro et al. 
2023). The combined limits are stronger than the results from 
the HESS observations of WLM dlir (Abdallah et al. 2021) at 
the DM masses above 10 TeV. Compared to other results 
given by analyzing dSphs, it can be seen that our constraints 
are weaker, mainly owing to the fact that dwarf irregular 
galaxies are more distant and thus have smaller J-factors. The 
advantage of our results, however, is that the J-factors of dIrr 
galaxies can be more accurately determined compared to 
dSphs and the constraints are thus more robust, as they are 
rotationally supported galaxies. In contrast, as pressure- 
supported systems, dSphs have larger systematic uncertain- 
ties in the J-factor determination (Battaglia et al. 2013). 
However, for dlrrs, there exists another systematic uncer- 
tainty from the determination of the boost factor, as discussed 
in Section 2. For this reason, here we also present the upper 
limits given by the joint analysis when the boost factor is not 
taken into account, and the results are shown as dashed black 
lines in Figures 3 and 4. 

From the individual constraints in Figure 3, we can see that 
the strength of the constraints depends mostly on the values of 
the dlrr J-factors, so the uncertainty of our results comes 
primarily from the J-factor uncertainty, which we discuss 
below. The (statistical) uncertainty on the dlrr J-factors 
includes mainly the error on the DM profile parameters 
obtained in the fit to the RC, which depends on the quality of 
kinematical measurements. This error is ~15% for po and ro 
and leads to a ~75% uncertainty in the astrophysical J-factor 
(Gammaldi et al. 2018). In addition to this, there are some 
systematic uncertainties in the J-factor estimation, e.g., the 
choice of the model for the DM density profile causes a 
systematic uncertainty of ~0.2 dex (Gammaldi et al. 2021). 
The J-factors used in this paper, which are derived in 
Gammaldi et al. (2018), assume a cored density distribution 
described by the Burkert profile (Burkert 1995), and such a 
profile is supported by the high-resolution HI observations (Oh 
et al. 2011, 2015) and can fit the rotation curves better 
(Gammaldi et al. 2021). However, DM-only N-body cosmo- 
logical simulations favor a cuspy Navarro—Frenk—White profile 
(Navarro et al. 1997). 

The largest uncertainty in the predicted DM annihilation flux 
comes from the uncertain nature of the dark matter substructure 
(i.e., the uncertainty in the enhancement factor B). As can be 
seen from Figure 4, whether the boost factor is considered or 
not, there is a factor of ~5 difference in the results. In the case 
of considering the boost factor, different choices of the 
parameters of the subhalo population (e.g., the slope of the 
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Figure 3. Upper limits on the DM annihilation cross section from the IceCube observations of dIrr galaxies for different channels. The black lines indicate the upper 
limits from the combined analysis of dIrrs. The two southern-hemisphere sources (WLM and NGC 6822) give weaker constraints especially in the lower mass regime. 


subhalo mass function, the minimal subhalo mass) lead to 
distinct (boosted) J-factors, as has been discussed in Section 2 
(see also Refs. Liang et al. 2017; Moliné et al. 2017; Gammaldi 
et al. 2021). Overall, for a dIrr of 10? Mo, the boost factor 
varies over a range of about 2.1—15.9 (Moliné et al. 2017), and 
thus the constraints (with B) on the DM parameters have the 
same level of uncertainty. 


5. Conclusions 


Pressure-supported dSphs are often considered as one of the 
most promising targets for indirect DM searches because of 
their close proximity and negligible astrophysical background. 
The problem, however, is that the DM content of dSphs and 
their astrophysical J-factors are usually subject to significant 
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Figure 4. A comparison of our constraints (black solid and dashed lines) with 
other existing limits in the literature. The results plotted for comparison are 
based on the observations from IceCube (Aartsen et al. 2013a), Fermi-LAT 
(Ackermann et al. 2015; Gammaldi et al. 2021), HAWC (Alfaro et al. 2023), 
HESS (Abdallah et al. 2021), MAGIC (Miener et al. 2022) and VERITAS 
(Acharyya et al. 2023). 


systematic errors due to their kinematic uncertainty. In contrast, 
dIrrs are rotationally supported dwarf galaxies, which are also 
DM-dominated objects. Unlike dSphs, the kinematics of dIrr 
galaxies provide a good estimate of their DM halo density 
distribution, and hence of their astrophysical J-factors. 
Although there exists certain star formation activity within 
dlrrs, it has been shown that the astrophysical background 
component is negligible (Gammaldi et al. 2018). They are 
therefore interesting novel targets for indirect detection of dark 
matter. In this work, we search for signals from UHDM 
annihilation in seven dIrrs using 10 yr of IceCube neutrino 
data. In the calculation of the J-factors of the dIrrs, we use a 
Burkert profile to characterize their DM distribution and 
consider the boost factors due to the existence of subhalos. 
Also, considering the angular extension of dlIrr galaxies in the 
sky, we analyze them as extended sources. 

We do not detect any significant signal of DM annihilation in 
the directions of these seven dIrrs. Based on the null results, we 
can place constraints on the model parameters of dark matter, 
since the emission from the annihilation should be detected if 
the annihilation cross section is large enough. Consequently, 
we derive the upper limits on the annihilation cross section at 
the 9596 confidence level for the DM with masses between 1 
TeV and 10 PeV. For the single-source analysis, the dlrr 
DDO154 gives the strictest exclusion limits. We also compute 
the results of the joint analysis of the five dIrrs in the northern 
sky. The joint analysis gives stronger (by a factor of 2-3) 
constraints on the annihilation cross section compared to the 
single-source analysis. 


Lii et al. 


In contrast to other experimental results, our constraints are 
equivalent to those given by the HAWC observations of dIrrs 
(Alfaro et al. 2023) and stronger than the constraints from the 
HESS observations of the WLM dirr (Abdallah et al. 2021). 
We also find that the exclusion limits of dIrrs are weaker than 
but close to those of dSphs (especially for the latest VERITAS 
results (Acharyya et al. 2023) at PeV masses), considering the 
smaller J-factor uncertainties of dlrrs, our results, therefore, 
provide a valuable complement to the existing results. 
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